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SYMBOLS 


Cauchy's  operator 
Differential  operator 
2 

D.D.  or  7 
1  x 

Strain  dyadic 
Body  force  field 
Galerkin's  vector 
Idemfactor 
Concentrated  load 
State  of  stress 
Stress  dyadic 
Surface  tractions 
Displacement  field 
Region  of  elastic  space 
Surface  of  V 
Convention  del  operator 

Determinant  of  C..  (D) 
lk 

Kroenecker  delta 
Permutation  symbol 
Lame/,s  constants 


Poisson's  ratio 


1.  Introduction 


The  equations  of  equilibrium  of  an  isotropic  linear  elastic 
body  V  are  called  Cauchy's  equations  and  may  be  written  as 

cik(D)  uk  +  PFi  =  0  tt-D 


where  C^^(D)  is  called  Cauchy's  operator  and  is  defined  below 

C.k(D)  =  p  D2  &.k  +  (\  +  u)  D.Dk  .  (1.2) 

The  Cauchy  equations  are  formulated  with  the  assumption  that  V  is  a  con¬ 
tinuous  medium.  This  ultimately  leads  to  the  restriction  that  u^  is  at 

least  piecewise  continuous  of  class  C^  throughout  V.  A  solution  to 
Eq.  (1.1)  which  satisfies  this  restriction  is  called  a  regular  solution*. 
The  form  of  a  regular  solution  to  Eq.  (1.1)  depends  upon  the  shape  of  the 
body  and  the  type  of  loading  imposed. 

Singular  solutions  exist,  however,  which  fail  to  exhibit  this  con¬ 
tinuity  at  a  point,  along  an  arc  or  over  a  surface  of  V.  The  purpose  of 
this  report  is  to  discuss  solutions  which  are  singular  in  some  deleted 
neighborhood  of  one  or  more  points  in  V.  Such  a  solution  has  one  or 
more  isolated  singularities  and  is  called  an  isolated  singular  solution. 

As  an  isolated  (nonremovable)  singular  point  x  in  V  is  approached, 
the  limits  of  u.  and  their  derivatives  do  not  exist.  In  fact,  the  dis- 

l 

placements  and  stresses  at  x  are  unbounded.  The  isolated  singularity 
may  be  considered  as  the  manifestation  of  some  type  of  load  or  moment. 

In  this  report  the  physical  significance  which  may  be  attached  to  various 
forms  of  isolated  singularities  will  be  discussed. 

It  will  be  readily  seen  that  at  least  three  useful  functions  are 
served  by  furthering  one's  understanding  of  isolated  singularities  of 
Eq.  (1.1).  There  are  numerous  cases  in  which  loading  is  closely  approxi¬ 
mated  by  a  concentrated  force  internal  to  a  body.  Also,  singular  solu¬ 
tions  may  be  used  as  influence  functions  to  obtain  regular  solutions. 

This  is  done,  for  example,  in  Betti's  [1|  adaption  of  the  method  of 
singularities  to  isotropic  elasticity.  Finally,  the  investigation  of 
point  defects  in  metallurgy  involves  these  solutions  (e.g.,  see  Simmons 
et  al  [2]). 


*Additionally ,  if  V  is  unbounded  it  is  required  that  the  u,  =  0(l/r) 

9  v 

and  u.  . =  0(l/r“)  uniformly  in  the  spherical  coordinates  -  ,  :  as  r 
J 

approaches  infinity. 
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First,  a  solution  of  Eq.  (1.1)  must  be  obtained.  Because  the 
Cauchy  equations  are  elliptic  (assuming  that  C..^(D)  *-s  positive  definite), 

techniques  analogous  to  those  found  in  potential  theory  may  be  used  f3'. 

In  Paragraph  2,  Cauchy's  equations  aia  written  in  biharmonic  form.  This 
is  possible  through  the  employment  of  the  Galerkin  vector  which  is 
developed  in  Paragraph  2. 


2.  Galerkin  Vector 

A  very  concise  derivation  of  the  Galerkin  vector  was  presented 
by  Somigliana  f 4 1 .  Choosing  a  transformation  K^(D)  so  that 


ui  *  Kii<<D>  vt<  ■ 


(2.1) 


the  Cauchy  equations  may  be  written  as 


Cik<D>  %m<D)  Vm  +  =Fi  '  0 


(2.2) 


Now  it  is  specified  that  Kj  is  the  matrix  of  cofactors  of  Cauchy's 
vector;  that  is 

K.  (D)  :  C.  C.  =  uD2  [*('.  +  2u)D2=  .  -  (\  +  n)D.D  1 

ip  2!  rjk  pmn  jm  kn  l  ip  l  pj 


By  definition 


(2.3) 


C..  (D)  K,  (D)  =  •  .  -(D) 

ikv  km  im 


(2.4) 


where  (D)  represents  the  determinant  of  C\^(D) 

(D)  =  ..2(  .  +  2n)  D6 


(2.5) 


Thus,  Cauchy's  equations  may  be  written  in  terms  of  the  variable  v., 
which  will  be  called  an  auxiliary  vector  function 


9  .  6 

“( .  +  2  :)  I)  v.  +  „*F.  =  0 
i  i 


(2.6) 
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Finally,  a  new  auxilliary  vector  function  is  defined  as 


8i  “7 


v. 

i 


(2.7) 


so  that 


^  =  U  ('■ 


2u)  D“ 


(k 


u)  D  D.g 
m  i  m 


(2.8) 


A  particular  solution  to  Eq.  (2.6)  may  be  found  by  using  the  concept 
of  influence  functions  which  are  analogous  to  Green's  functions  in  poten¬ 
tial  theory.  An  insight  into  this  approach  can  be  gained  from  the  one 
dimensional  string  problem.  Consider  a  string  L  units  in  length  which  is 
fixed  at  both  ends  and  which  lies  along  the  x^  axis.  If  the  string, 

under  a  tension  t,  is  subjected  to  a  transverse  unit  force  applied  at 
x^  =  s,  the  equation  of  equilibrium  is 


tx 


2,  11 


(2.9) 


f>(x^  -  s)  is  the  Dirac  delta  function  defined  by  Chen  [5] 


|  0,  for  x^  •/.  s 
|  <v ,  for  x^  -  s 


A  solution  of  Eq .  (2.9)  indicating  the  vertical  displacement  of  the 
string  is 


•" 2  K(xr  s)  (2.10) 

where  K(x^,  s)  is  called  an  influence  coefficient.  Importantly,  Eq.  (2.10) 

can  be  used  to  solve  problems  concerning  distributed  transverse  loadings 
q(x^)  because  it  can  !x-  shown  thac  the  resulting  equilibrium  equation 

has  the  solution 

L 

-  J  K(Xp  s)  q(sl  u«  .  <2.  Ill 

t'» 
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A  solution  to  Eq.  (2.6)  which  corresponds  to  a  body  force  pF^  can 

be  determined  in  an  identical  manner.  To  avoid  complications  introduced 
by  finite  regions,  a  body  V  of  infinite  proportions  is  chosen.  Imagine 
a  body  force  pF^  which  is  only  nonzero  over  a  finite  region  V.  Then, 

if  S(x,  |)  is  the  solution  of 

A(D)  S(x,  i)  +  S3(x  -  g)  =  0  (2.12) 

where  “  0  is  the  volume  Dirac  delta  function,  then  a  solution  to 

Eq.  (2.6)  at  a  point  x  is 


O  pF.(5)  dD(|) 


(2.13) 


Fritz  John  [6]  presents  the  general  form  of  S(x,  fj).  The  particular 
form  of  S  (x.,  f- )  for  Cauchy's  equations  was  found  to  be  [7] 


S(x,  jp  =aR(x,  P) 


(2.14) 


where 


a  = 


4 *n2(\ 


+  2p) 


R(x,  0  = 


(2.15) 


Thus,  a  particular  solution  to  Eq.  (2.6)  is 


k=a  f' PFi(l 


(4)  R(x,  O  dD(4) 


(2.16) 


In  this  form  the  Galerkin  vector  q^  in  Eq.  (2.8)  yields  a  solution  u^ 

which  represents  che  displacement  within  an  infinite  linearly  elastic 
body  V  because  of  a  body  force  field  pF.  acting  only  within  a  finite 
region  V.  L 
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3. 


Limit  Solution 


The  purpose  of  this  section  is  to  employ  the  limit  solution 
technique  of  Sternberg  and  Eubanks  [8]  to  develop  and  discuss  various 
types  of  isolated  singularities  in  linear  isotropic  elasticity. 

-Again,  entire  three  dimensional  space  \T  is  considered.  This  time, 

(n) 

imagine  a  sequence  of  regions  {v^  }  each  containing  the  origin  and  a 

corresponding  sequence  of  body  forces  fpF.^)  such  that  pF.  ^  acts 

/  \  y  ^  \  /  ^  \  X  X  /  ^  \ 

over  region  V'  ,  pF.  'over  v  ,  etc.  It  is  required  that  each  pF.  'be  of 

class  C  in  V°°  and  be  zero  at  any  point  outside  of  V'  .  It  has  been 

shown  that,  under  these  circums  ances,  a  particular  solution  u.^(x) 

/  \  /  \  1 

exists  for  each  V'  '  and  corr  ponding  body  force  field  pF^  .  In 
addition,  as  n  approaches  infinity  the  length  of  the  longest  chord 

/  n  \ 

within  V'  '  must  uniformly  approach  zero. 

The  transformation  represented  by  Eq.  (2.14)  is  now  employed, 
yielding 


,(n) 


(3.1) 


Thus,  if  f.  is  a  typical  point  within  V'  1  and  x  is  any  other  point  within 
V‘  , 


;.(n)(x)  =  o  Jr(x,  «)  pF.(n)(0  dv(n)(0 


(3.2) 


The  limiting  case  of  Eq .  (3.2)  is  defined  as 


^(x)  =  Lim  g.(n)(x)  =  Lim  a  /  R(x,  ?)  pF.  (n)  (s)  dV(n)(  = 
n  ->cv  n  — »>  j 


)  •  (3.3) 


Richardson  f 7 ]  showed  that,  in  general, 


(x)  ,  .  (n) 

G.  ..  =  Lim  g.  . 
L» dk  n  .  „  i.jk 


(3.4) 
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Therefore,  U^.  computed  from  G,  represents  the  limiting  solution,  a  dis¬ 
placement  field  with  an  isolated  singularity  at  the  origin.  The  forms 
which  this  solution  assumes  will  be  discussed  in  later  sections. 

R(x,  ■)  may  be  expanded  in  the  Taylor's  series  for  three  variables 
shown  below 


R(i,  0=  £ 

i ,  j , k=0 


l-..1.)1+j+k  DiDJDk  (r)  -iJ 

i!j!k!  ulD2°3  1  23 


r  =  r(x) 


(3.5) 


Thus,  Eq .  (3.3)  becomes 


G  (x)  =  *. 


E 

i ,  j  ,  k=0 


.  i  \  i+ j+k  .  .  , 
i! j!k!  D1D2D3 


Lim  /  oF 
n--.  •  J 


F(n)0) 


.  i  J  „  k 


]  2  3 


dV(n)(’)  , 

(3.6) 


Note  that  the  summation  sign  may  bo  removed  from  under  the  integral  sign 
because  Taylor's  series  is  uniformly  convergent.  Defining  the  moment 

integrals  which  are  of  order  i-i-j+k  and  operators  '  by 


MiJk  -  I, im  Lr 
n  . .  J 


(nV)  uv(nV) 


(3.7) 


.  ijk 


. . i+j+k  .  .  , 

- -  D^D^D. 

i!j!k!  1  -  J 


(3.8) 


where  i+j+k  =  p  which  is  the  order  of  the  operator  1  allows  one  to 
write  the  Galerkin  vector  as 


0< 


i  jk 


« r ) 


'  j,U(> 


( 3 .  J  i 


Thi°  is  t  lorle.  hick  vnivv'ts  tin.  \»I.  ■*.; 

"  ar  ties  et  >l|  n-doi\>.  Add  it  \>nal  !'•.  by  cho 
c;ii  cries.  >ne  o  oi  -.rile  the  ..aleck  in  vector 
1  oi":v;  ef  i. m!.uc:  singularities. 


•  vc»  t"i‘  "or  iso ■>: nuu- 
■  apovon'-iatc  con.,  of 
CO'*r  '’Htl'n  i  !'•  t'  sivt  i  lie 


x  i  lc 

It  is  of  interest  to  investigate  the  nature  of  M  .  Note  that 
if  the  sequence  of  regions  and  the  corresponding  body  force  fields 

pF.  are  symmetric  about  all  axes,  then  the  integrals  with  an 

^  ( n ) 

odd  i,  j,  or  k  are  zero.  The  shape  of  the  region  V  is  arbiti*ary  so 

symmetry  can  be  assumed;  however,  pF.^n')  are  not  generally  symmetric. 

Later  in  this  report  some  manifestations  of  symmetry  of  pF^  '  will 
be  discussed. 

In  Paragraphs  4,  5,  and  6  it  will  be  shown  that  Eq.  (3.S)  may  be 
written  in  the  form 


G.  (x)  =  o  (r)m.  +  6  .  (r)m.  (r)m.  -!■... 

i  -  i  j  ij  jkv  ijk 


(3.10) 


4.  Fundamental  or  First  Order  Singularity 

The  term  in  Eq.  (3.10)  corresponding  to  i=j.-k.  0  is 


Gt(x)  -or 


Lim  f pF.(nV 
n_»  oo  J 


)  dV(n)(f-) 


(4.1) 


From  Eq.  (2.8)  it  is  seen  that  the  displacement  vector  is  expressed  in 
terms  of  the  second  derivatives  of  g^.  Therefore,  if  g^  is  of  0(r)  as 

r  ->  0  then  u^  will  be  of  0(l/r)  r  _>  0.  This  case  shall  be  called  the 

first  order  singularity.  Superscripts  on  G.  will  indicate  the  order  of 
singularity.  Defining 


(r)  -or,  m.  -  Lim  /pF.(n)0)  dV("V) 
n  sc  J 


(4.2) 


produces  the  first  term  of  Eq.  (3.10) 


G:(x)  -  <  (r)m. 


(4.3) 
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It  is  readily  noted  that  Eq.  (4.2)  may  be  interpreted  as  the  Galerkin 
vector  leading  to  a  solution  which  corresponds  to  a  concentrated  load 
acting  at  r  =  0  if  it  is  specified  that 


Lim 

n  ->eo 


/ 


pF.(n)(0  dV(n)(|) 


L. 

1 


(4.4) 


In  the  following,  the  fundamental  displacements,  displacement  gradients, 
strains,  and  rotations  are  given.  Throughout  the  remainder  of  this 
section,  i?;j*k  and  i,  j,k  will  not  be  summed. 

The  strain-displacement  relations  are  given  by 


U1  = 
i 


—  - y  4(1  -  v)  m.x3  +  (3  -  4v)  m.x3 

-  L  11  1  J 


8rtu(l  -  v)r“ 

+  (3  -  4v)  nux3  +  nux^  +  nyc.XjJ  •  (4.5) 

The  displacement  gradients  are 

..1  1  I",  ,  3  3  3  „  2  2 

b.  .  -  - s-  4(v  -  1)  m.x.  •!•  m.x.  m,  x,  -  2m.x.x.  -  2m,x.x. 

1'1  8«u(l-rtr5  L  11  ”  'l 

2  2  2  ?  1 
+  -  1)  m^Xj  +  (4v  -  1)  mixixk  +  siux^  +  •  (4.6) 

However,  for  the  off-diagonal  terms, 

u}  -s  - - - r  [m.x3  +  (4v  -  3)  m.x^  +  (4v  -  6)  m.x3x.  -  2m.x.x3 

l>-  8  „U  -  v)r5  L  1  1  1  J  1  1  J  J  1  J 


2  2  1 
t  m.x.x,  -i  (4v  -  3)  m.x.x,  -  3m,x.x.x. 

J  i  k  i.]k  k  i  j  kj 


Next,  the  strains  are  found  to  be 


(4.7) 


e!. 


if  3  3  3  o  2 

- c  4(i  -  1)  m.x.  -im.x.  .  nix,  -  2m.xTx.  -  2m,x.x, 

L  ii  JJ  k  k  j  i  j  k  i  k 


8"u(l  -  v)r 


(4'.  -  l)  m.fx.x3  +  x.x,3)  •'  m.x.x,4"  m,x3x,l 
i\  1  J  i  k/  j  j  k  k  j  kj 


(4.8) 
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For  the  off-diagonal  terms, 


E1. 


1J  8>cix(l  -  v)r 


- - -  (2v  -  1)  m.x^  +  (2v  -  1)  m.x^  +  (2v  -  4)  nux?x. 

—  I/'t  T~ 


J  1 


i  J 


i  i  J 


+  (2v  -  4)  m.x.xf  +  (2v 
J  i  J 


9  2  11 
1)  m.x.x.“  +  (2v  -  1)  m.x.x,  -  3m,x.x.x,  =E.. 

jik  i]  k  k  l  j  kj  ji 

(4.9) 


Finally,  the  rotations  are 


Ol  r 


1J  8rt|i(l  -  v)r 


- - =■  4(1  -  v)  m.x?  +  4(v  -  1)  dux'?  +  4(i  -  I)  m.x?x. 

.  ,a  5  v  J  i  ij  1  1  J 


2  2  21  1 
+  4(1  -  v)  m.x.x.  +  4(1  -  v)  m.x.x,  +  4(v  -  1)  m.x.x,  =  -Q.. 

J  i  J  v  J  i  k  '  i  j  kj  ''ji 


(4.10) 


So,  the  displacements,  displacement  gradients,  strains,  and  rota¬ 
tions  for  a  state  of  stress  with  a  first  order  isolated  singularity  at 
the  origin  have  been  defined  in  terms  of  m.,  it,  v,  and  the  coordinates 
of  thf  point. 


a  i  ""r 

A.  .  v 
i-J- 


It  •>’  now  possible  to  obtain  higher  order  isolated  singularities  by 
.i(ir  superposition  technique.  Flexibility  influence  coefficients 
£)  can  be  used  to  describe  the  effects  of  a  force  Q./V  acting  at 


a  po.  t  f  on  the  displacement  u.  at  some  other  point  x  (Figure  1).  In 


fact,  the  displacement  is  simply 


ujoo 


(4.11) 


Under  conditions  shown  in  Figure  1,  Eq.  (4.11)  will  be  identical  to 
Eq.  (4.5).  This  represents  the  fundamental  or  first  order  singular 
solution.  Similarly,  the  displacement  at  x  because  of  a  concentrated 
force  -  Q./c  acting  at  f  +  e  is 


A .  . (x ,  r 
ij  “  - 


c)  ^ 

—  c 


(4.12) 
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Because  the  force-displacement  relationship  is  considered  linear,  the 

displacements  created  by  Q./c  at  and  -Q./c  at  ■"  c  are  simply  the 

J  r 

superposition  of  Kqs.  (4.11)  and  (4.12)  (Figure  2).  Taking  the  limit 
of  the  resulting  expression  yields  a  second  order  singular  solution 


V5,  :>  - 


9  f  1  1  1  Ai 

U7(x)  =  Lim  U.  (x)  +  U.  (x)  -  him  — 

.%,0  [  1  1  I  r.,0 


which  is  recognized  as  the  first  derivative  of  U.(x), 


(4.13) 


Specifically,  in  Figure  3,  Q  is  parallel  to  the  x^  axis  and  c  is 
measured  along  the  x^  axis.  Therefore,  let  u!"(x)j  and  u!“(x)|  be 

represented  by  Eq.  (4.5)  where  only  is  nonzcroT  For  convenience,  the 

point  ■’  is  designated  at  the  origin.  Because  the  influence  coefficients 
are  dependent  upon  the  v>'ct,v»*  R«x,  '),  the  displacements  at  x  caused  by 
the  forces  shown  are 


10 


x3 


Figure  4.  Force  Couplet 

*Three  force  doublets  of  equal  magnitude  acting  along  the  x^,  x?,  and 
x,  axes  comprise  a  center  of  dilatation. 
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The  generalization  of  this  process  is  intuitively  simple.  Doublets 
and  couplets  of  any  magnitude  may  be  generated  in  various  directions  or 
may  themselves  be  superposed  to  form  third  or  higher  order  singularities. 
This  procedure  is,  however,  artificial  and  tedious.  In  Paragraphs  5  and 
6,  the  limit  solution  technique  will  be  employed  to  form  expressions 
for  all  possible  second  and  third  order  isolated  singular  solutions. 


Second  Order  Singularities 


100  010 


It  is  now  observed  that  the  terms  >I>(r),  ‘J>(r),  'l'(r)  and  M 

0,  M^^  can  be  treated  as  components  of  first  and  second  order 
tensors.  Making  the  following  definitions 


->  (r)  =  -  ,  m  =  Lim  f  pF  (n)(|)  P  dV(n)(*)  .  (5.1) 

J  -1  n  ->  <xJ  i  - 

The  Galerkin  vector  corresponding  to  second  order  singularities  may  be 
written  as 

G?(x)  =  <t>.(r)  m..  (5.2) 

i-  j  ij 

which  is  identical  to  the  second  term  of  Eq.  (3.10). 

To  examine  the  physical  characteristics  of  various  second  order 
singularities,  one  may  first  allow  all  nu.  except  m^  to  be  zero. 

In  this  case,  Eq.  (5.2)  becomes  1‘1 

G2(x)  -  Jo^(r)  mn,  0,  0]  .  (5.3) 

One  may  say  that  the  resulting  displacement  pertains  to  a  force  doublet 
along  the  x^  axis.*  Similarly  if  ir^  0 r  m^  are  the  only  nonzero  com¬ 
ponents  of  m. j ,  then  the  doublets  are  along  the  x^  or  x^  axis,  respec¬ 
tively.  The  significance  of  the  off-diagonal  terms  is  shown  by  setting 
all  nKj  equal  to  zero  except  •  Now  the  Galerkin  vector  is 


*If  only  m^,  m^,  are  nonzero  and  m^.,  -  -  nu^,  a  center  of 

dilatation  is  produced. 
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G“(x)  =  [o,  4-3(r)  m23>  o] 


(5.4) 


This  coi-responds  to  a  force  couplet  about  the  axis  with  the  forces 
•  .1  ♦ _ * _ . 


in  the  x2  direction. 

If  pF.  ^  are  symmetric  about  all  axes  then  all  m„  are  zero  and, 
hence,  there  are  no  second  order  singular  terms.  If  pF_.(n)  are  sym¬ 
metric  about  one  or  two  axes  then  certain  components  of  the  tensor  m. . 

( n)  ^ 

will  be  zero.  For  example,  if  pF.  is  symmetric  about  the  plane  per¬ 
pendicular  to  the  x^  axis  then  =  m-3  =  0  and  the  resulting 

Galerkin  vector  is 

G2(x)  ={[«’1(r)m11  +  4s 2 ^mi 2] 5  [<>i(r)m2i  +  *2(r)m22]»  [<1>l(r',m3l +<’2^r')m32]} 

(5.5) 

This  indicates  doublets  on  the  x^  and  axes  as  well  as  couplets  about 
the  x^,  x2,  and  x3< 

If  pF^n^  are  symmetric  relative  to  planes  along  two  axes  (for 
example,  the  x9  and  x3  axes)  then  only  the  m^,  n^,  and  m3^  components 
are  nonzero.  Hence,  the  Galerkin  vector  becomes 


G2(x)  JV1(r)m11,  <t1(r)m71,  -.^(rjr.^j 


(5.6) 


which  corresponds  to  a  doublet  on  the  x^  axis  and  couplets  about  the 
x9  and  x3  axes. 

Every  second  order  tensor  may  be  expressed  as  the  sum  of  a  symmetric 
and  antisymmetric  tensor.  Therefore,  may  be  written  as 


S  A 
m.  .  m.  .  -i  m.  . 

ij  U  ij 


where 


(5.7) 


SI/  \  S 

11.  .  -Ini.  .  -i-  m. .  ,  m,  . 
T-J  2\  ij  J x/’  IJ 


(5.8) 


A  1  /  \  A 

1.  .  -Im.  .  -  m. .  ),  m.  . 
ij  2\  ij  31/ *  J.j 


(5.9) 
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Now,  the  products  of  <K(r)  ar»d  are  given  below: 


„2S.  .  .  .  S  a  S 

G.  (x)  -  $ . (r)  m. .  - - x.m.. 

i  —  j  ij  r  j  xj 


(5.10) 


Similarly,  the  products  of  <^(r)  and  m£j  are 


~2A  ,  .  A  /  s  A  a  A 

G.  (x)  =  4>  .  (r)  m..  =  -  -  x.m.. 

i  -  J  ij  r  j  ij 


(5.11) 


The  displacement  fields  corresponding  to  Eqs.  (5.10)  and  (5.11)  are 

„2S  1  ( T,  ,  ,,  S  S  Si  3  „  _  S  3 

U  =  - 5  4(y  -  1)  m  .  +  m.  +  m..  x  -  (4v  -  2)  m  x 

16jr|i(l  -  v)r  (L  11  JJ  kkJ  1  1J  J 


//  ON  S  3  .  .  ..  S  2  ..  S  2 

+  (4v  -  2)  m..x.  +  4(r  -  1)  m.  .x.x.  +  4(\  -  1)  m.,x,x, 

lk  k  ij  i  j  lk  i  k 


'  /  S  .  S  ,  S  1  2  [*..  ..  S 

+  (4v  -  1)  m..  -  2m..  +  m.  ,  x.x.  4  (4y  -  1)  m.. 

xi  JJ  KkJ  x  J  v  1.1 


S  S 

. .  4m..  - 

ii  JJ 


S  1  2 

kk  l  k 


+  <4v  - 2)  4xK +  <4v  - 2)  muxfk  -  6"lkxixjxk 


—  no  sum,  i  /  j  ^  k 


(5.12) 


_  - - - —  |  4(1  -  . )  nu.x?  •*  4(1  -  »  J  nAx?  i  4(1  -  )  nu.xTx. 

1  16i»u(l  -  :)r5l  1Jj  lkk 


4  4(1  -  \)  nn.  x?x.  «  4(1  -  i)  nA  x?x.  ^  4(1 
lk  i  k  v  lk  j  k 


)  m.  .x  ,x“ 
ij  J  k 


- no  sum,  i  -/  j  /  k 


(5.13) 


It  is  easily  shown  from  Eqs.  (4.8)  and  (5.12)  that  if  m.  is  the  only 

Q 

nonzero  component  in  Eq.  (4.8)  and  if  m  =  m  and  is  the  only  nonzero 
S  ii  i 

element  in  the  m..  tensor,  then 
ij 


2S  1 

U1  E11 


(5.14) 
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This  can,  of  course,  be  generalized  for  the  second  and  third  components 

2S  2S  1 

of  U  as  well.  Thu;.,  under  these  conditions  U.  and.  E..  correspond  to 

a  force  doublet  in  the  x.  direction.  1 

1 


S  S 

m^  in  Eq.  (4.8)  and  =  m ^ 

and  are  the  only  nonzero  terms  in  nu  . ,  then 

ij 


Additionally,  if  =  m2 


S 

m33 


(5.15) 


and  represent  a  center  of  dilatation.  Farther,  from  Eqs.  (5.10)  and 
(5.13)  it  can  be  shown  that  each  component  of  the  antisymmetric  part  of 
the  second  order  singularity  is  identical  to  the  rotation  vector  about 

A 

its  axis.  That  is  (if  the  signs  of  m. .  are  dictated  by  the  right  hand 
rule),  1J 


i  *  j  •/  k 


and  corresponds  to  a  force  couplet  about  the  x^  axis. 


(5.16) 


6.  Third  Order  Singularities 

The  third  term  in  Eq.  (3.10)  can  also  be  expanded  and  special 
ized  for  particular  cases.  The  expansion  is  shown  as 


Ci<S>  =Vr)  mijk 


(6-1) 


where 


m  ~  Lim  /pF  (n)(0  dV(n)<0.  (6.2) 

J  n-.  co  J  J 


Only  one  particular  form  of  third  order  singularity  will  be  presented; 
however,  it  is  obvious  that  one  can  easily  express  any  specific  type. 
For  example,  the  following  Galerkin  vector  corresponds  to  a  singularity 
at  the  origin  which  can  be  called  a  double  center  of  dilatation  on  the 
x^  axis  if 


Vr)  0 


1  .  .  x  .x 

r  3 
r 


16 


(6.3) 


G  (x)  _  [<J,13(r)m131>  * 23 ^r)m232’  <>33^r^m333] 


where  2m131  =  2m232  =  m333. 


7.  Spherical  Biharmonics 

The  series  represented  by  Eq.  (3.10)  provides  the  foundation 
for  a  development  analogous  to  that  of  the  class  of  Newtonian  potentials 
called  spherical  harmonics.  It  may  be  recalled  that  a  Newtonian  p.  '-.en- 
tial  is  a  solution  to  Laplace's  equation  at  points  external  to  the  body. 
For  a  body  of  density  a (§  )  occupying  the  region  V  +  in  space,  the 

Newtonian  potential  at  x  is 


N<5>  "/ R0T) 


dV(" ) 


(7.1) 


By  expanding  1/R(x,  f)  about  x  in  a  Taylor  Series,  N  may  be  written  as 


E  ^(t)^ 


jk 


i»  j,k:-0 


where 


(7.2) 


Bijk  (?)  -  nw1  D}DiD5B) 


(7.3) 


■  <7-« 

It  is  easily  seen  that  B^j^(l/r)  is  homogeneous  of  degree  - (i+j+k+1) . 

Further,  1/r  is  harmonic.  Therefore,  each  B.  ..  (1/r)  must  be  harmonic 
2  LJk 

because  ' ’  operates  only  on  x^.  Thus,  by  definition,  (1/r)  are 

spherical  harmonics  of  degree  -(i+j+k+1). 
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A  comparison  of  Eqs.  (7.2)  and  (3.9)  enables  one  to  make  several 
observations.  First,  Eq.  (7.2)  represents  the  potential  function  of  an 
attractive  force  on  x  by  the  mass  occupying  region  V,  while  Eq.  (3.9)  is 
a  potential  function  of  the  displacement  at  a  point,  x  resulting  from  a 
force  or  moment  acting  at  the  origin.  It  was  shown  chat  the  coefficient 
Bij^(l/r)  are  spherical  harmonics  and  by  similar  reasoning  it  is  noted 

that  (r)  are  biharmonic  and  homogeneous  of  degree  (i-i-j-k),  a 

function  which  will  be  called  a  spherical  biharmonic  of  degree  (i-i-j-k) 


Thus,  with  the  analogy  clearly  established,  an  interesting  point 
may  be  investigated  by  first  noting  that  the  expansion  of  Eq.  (7.1) 
yields 


+  0 


+ - 


(7.5) 


All  terms  except  the  fundamental  singularity  may  be  eliminated  by 
prescribing  that 


i  .J,k 

12  3 


dV(') 


=  0 


for 


i  ■'  0,  j  '  0,  k  /  0 


(7.6) 


This  step  is  sufficient  to  insure  that  N  is  the  fundamental  singularity, 
but  is  it  necessary?  One  suspects  that  it  is  not  when  considering  the 
following:  A  complete,  homogeneous,  polynomial  of  degree  n  contains 
(n  +  l)(n  +  2) / 2  coefficients,  only  2n  -F  1  of  which  are  arbitrary  f 9 ) 

? 

(assuming  that  all  coefficients  are  harmonic).  For  example,  the  0(-7) 
terms  are 
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^5  (&1  "  4 ' 4)  fa  1  dn0  +  (“xl  +  2x2  ■  4)  fa 2  dV(D 

+  ^-x2  -  *2  +  2x^)  y crf=3  dV(p)  +  6xjX2  J o?^  dV(0 
+  6x^2  JctlS3  dV^lN  +  6x2X3  J*P2*3  dv j  *  1 


(7.7) 


The  second  degree  complete,  homogeneous,  harmonic  polynomial  inside  the 
brackets  has  six  terms.*  The  fact  that  there  are  dependent  coeffi¬ 
cients  in  terms  of  order  greater  than  1  suggests  that  conditions  other 
than  Eq.  (7.6)  will  result  in  the  elimination  of  higher  order  terms. 
This  lack  of  uniqueness  can  be  illustrated  by  allowing 


Jo  =  2  dV(_)  =  fa*  d\’C)  =  fa 2  dvn 


(7.3) 


in  Eq.  (7.7).  In  this  case  the  first  three  terms  in  Eq.  (7.5)  are  zero, 
2 

The  0(0  terms  may  bo  written  in  powers  of  x^;  that  is, 

—j  (2Ml  -  M2  -  M3)  Xj  +  (“M1  +  iM2  ”  ^3)  x2  +  ("Ml  "  M2  +  "’‘'b)  x3 


where 


MpXjX,  +  M13XjX3  -)•  M?3x9x3 

M..  -  d\'  ( ’ ) , 

U  J  O  - 


(7.9) 


1  '  J 


(7.10) 


‘The  bracket ted  expression  may  he  discussed  with  generality  because 

'’n  +  1 

if  B  is  a  spherical  harmonic,  then  B  /r“  is  also  a  spherical 
harmonic  f9). 
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and 


M. 

1 


/■ 


0f-t.  dV(|) 


Now  the  first  three  terms  of  Eq. 


(7.9)  are  zero  if 


'  2  -1  -1' 

fM}[K]  - 

^2 

-1  1  -1 

r3 

-1  -1  2_ 

(7.11) 


(7.12) 


Thus,  a  condition  other  than  Eq.  (7.6)  exists  which  eliminates  the 
2 

0(?.)  terms. 

l 

3 

The  0(>‘.)  term  has  10  coefficients.  One  can  write 
x 


+ 


_/-(-3’ i  -  }‘A+ UfA)  dv<a>  v32+  /  H  -  A- 1  - 3  y;s)xv<±>  ^ 

fo  (-3’1  +  12=2. -  3  V2)  d»(0  m2x2+  J- (-3- j  -  3  Y 2+  12  Y 3) dv<.I>  *2*3 
J„(2-32-3Vl2-3V2)<IV(1)x3+  Ja(-  3-3+  3^2  -  3  V2)  dv(.)  x2x. 

j"{-A-3’A  +  n‘A)'1V'tJ*h  +  J  30~  ■ 


(7.13) 


20 


It  may  be  seen  that  the  entire  term  vanishes  if  all  integrals  are  zero, 
or  if 


Jf  fi  dv^)  =  /Vj  dv(l>  =  jrt x  k  dV(-}  =  ai>  *k 


{^1*2*3  dV(=)  =  0 


(7.14) 


and  so  on  for  all  orders. 

To  visualize  the  theoretical  possibility  of  such  an  occurrence, 
imagine  a  sphere  of  radius  c  centered  at  0  and  a  spherical  shell  with 
an  outer  radius  a  and  an  inner  radius  b  also  centered  at  0.  If  b  c 
the  sohere  is  encapsulated  inside  the  shell.  The  potential  of  attraction 
of  the  sphere  relative  to  a  point  x  lying  outside  the  sphere  is 


c 

so -  7^7  /  dr<:>  -  5^i>  '3  ■  <7-15> 

0 

Similarly,  the  potential  of  the  shell  relative  to  a  point  x  lying  inside 
of  the  inner  wall  is 

a 

Nt  =  '  J  r(f)  dr(-)  =  2^' (a2  -  b2)  .  (7.16) 

b 


Of  course,  the  potential  of  a  point  lying  between  the  sphere  and  shell 
is  N  »  Nq  +  Nj . 

If  one  prescribes 


.  i  _ 


(7.17) 


then  N  is  zero  without  requiring  that  all  integrals  vanish. 


’  ‘  ’  7'TV’®I7!  rT^V f ’ 


The  preceding  investigation  was  conducted  with  a  harmonic  function 
f.  Its  results,  however,  can  be  applied  to  a  biharmonic  function  g 
by  noting  that  an  arbitrary  harmonic  function  8  can  be  defined  so  that 


^2g  =  n 


(7.18) 


Thus,  in  the  expansion  of  the  biharmonic  vector  there  will  be  depend¬ 
encies  among  coefficients  beginning  with  the  0(’l\)  term,  because  this 

2  1 
corresponds  to  the  0(*^)  term  in  the  expansion  of  . 


8.  Betti's  Reciprocal  Work  Theorem 

In  Paragraph  2,  Galerkin's  vector  £(x^)  was  used  in  discussing 

the  development  of  a  particular  solution  of  the  Cauchy's  equations.  In 
general,  these  functions  are  analogous  to  potential  functions  and,  as 
stated  previously,  several  techniques  and  theorems  of  potential  theory 
are  applicable  in  the  treatment  of  the  equations  of  elasticity.  A 
lucid  explanation  of  the  application  of  potential  theory  to  the  solution 
of  elliptical  equations  is  found  in  Courant  and  Hilbert  [3] .  One  of 
the  most  important  aspects  of  potential  theory  is  the  theorem  and  iden¬ 
tities  of  Green.  Betti  fl],  in  1872,  adapted  Green's  identities  to 
elasticity  and  thus  formulated  the  reciprocal  work  theorm.  In  this 
section  the  reciprocal  work  theorem  and  a  particular  application  will 
be  discussed. 


Consider  two  stress  states  S(x)  and  S'(x)  which  are,  for  the  moment, 
assumed  to  be  regular  throughout  an  elastic  body  V.  If  the  unprimed 
state  has  a  zero  body  force,  the  equations  of  equilibrium  For  S(x)  and 
S '  (x )  arc 


h  ...  u  ,  .  =  0 
i.l kl  k,  Ij 


(8.1) 


H. .. .  u,'  . .  +  *r:  =  o 

i;kl  k,lj  i 


(8.2) 


Multiplying  Eq.  (8.1)  by  u|  and  Eq.  (8.2)  by  u^,  respectively,  and  sub¬ 
tracting  the  resulting  expression  yields  a  relationship  between  the 
primed  and  unprimed  states.  Integrating  this  relationship  over  V  results 
in  the  following: 


III.  ..  .  iu.  ,  .  u!  -  u,'  ,  .  u.  (IV  frl'lu.  dV  = 
I  ljkl  I  k,lj  i  k, ij  i  /  ii 


(8.3) 


Finally,  employing  Green's  theorem  and  noting  that 


M 


n . 
J 


u!  d-V 
1 


T .u!  d"V 

J  1  1 


one  may  write  the  reciprocal  work  theorem  as 


(8. A) 


dV  =  0 


(8.5) 


Imagine  now  that  lying  within  the  region  V  is  the  sequence  of 
regions  f V ^ n ]  having  the  properties  described  in  Paragraph  3.  Suppose 
further  that  the  primed  stress  state  corresponds  to  S^n  (x),  where 
pF.  (n\x)  is  zero  for  x  not  in  For  the  region  V^n\  Bq.  (8.5) 

becomes 


j[ T.u.<n)  -  T.<"V]  d  -V00  - 


dV(n)  =0  .  (8.6) 


Similarly,  applying  Bq.  (8.5)  over  the  region  V  -  yields 

-  J['riu.(n)  -  T.(n\i.  j  d  -V  +  J  Fiu.(n)  -  T.(nVj  d  V(n) 


+ 

Note  that  the  tractions  over  V ^ 11 '  as  well  as  must  be  considered. 

The  terms  T . ( 11  and  T.  involve  an  outer  normal  in  Bq.  (8.6)  and  an  inner 

normal  in  the  second  integral  of  Bq.  tS.7);  there  tore,  the  signs  are 
opposite,  fhe  addition  of  these  two  equations  results  in 


(n> 


u.d 

x 


('■- 


,Cn)\  _ 


=  0 


(8.7) 


0 


( 8 . 8  > 


Expanding  ik  about  the  origin  in  a  Taylor's  series  and  taking  the  limit 
as  n  approaches  infinity  yields 


0  =  Lim  f  fau^1^  -  T^Vl  d^v  +  ui^°)  Lim  f  PFi(n'>  dv^n) 
n  — >co  J  L  J  n  —>co  J 

+  u.  .(0)  Lim  /  oF.(n)f,  .dV(n)  +  u.  .,(0)  Lim  /  pF.  .F,  dV(n) 
I  i  i  l, ikv  /  '  i  j  k 

’ J  n  ->«  J  J  J  n~>K  J  J 


+  — 


(8.9) 


or 


0  -  J  [t.U,  -  TjU.] 


dAV  +  u.  (0)m.  +  u.  .  (0)m.  .  +  u.  (0)m.  ..  +  — 
ii  i,J  ij  i, ]kv  ljk 


(8.10) 


where 


r 


i 


Lim  T. 

l 

n->«- 


(n) 


And  so,  by  choosing  the  appropriate  order  singularity  (by  choosing  values 
of  nu,  m„ ,  etc.),  the  displacements  or  any  order  derivatives  of  the 

displacement  field  of  a  regular  stress  state  acting  throughout  can  be 
reproduced  if  the  displacements  and  tractions  are  known  on  "'V.  Note 
that  by  choosing  a  particular  isolated  singular  solution  U.,  one  corre¬ 
spondingly  chooses  the  form  of  nu,  m^ . ,  m^.^,  etc. 


9.  Discussion 

As  stated  previously,  the  object  of  this  report  is  to  offer  a 
physical  interpretation  of  various  first  and  higher  order  isolated  singu¬ 
lar  solutions  of  Cauchy’s  equations  using  Galerkin's  vector.  These  singu¬ 
larities  are  constructed  by  means  of  a  limit  solution  proposed  hy 
Sternberg  and  Eubanks  (8],  Some  forms  of  first,  second,  and  third  order 
isolated  singularities  are  investigated  and  will  be  followed  by  a  discussion  of 
uniqueness  and  an  application  of  Betti's  reciprocal  theorem.  In  the 
following  paragraphs,  various  salient  points  will  be  reviewed. 
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In  Paragraph  2  Cauchy's  equations  are  rewritten  in  a  form  more 
easily  solved.  This  form  is  a  nonhomogeneous  biharmonic  equation  in 
terms  of  an  auxiliary  vector  called  Galerkin's  vector.  The  solution  is 
derived  by.  employing  a  technique  analogous  to  one  used  in  potential 
theory. 

Next,  the  limit  solution  is  discussed  and  a  series  is  developed 
which  contains  all  first  order  isolated  singular  solutions  to  Cauchy's 
equations  in  the  first  set  of  terms,  all  second  order  singularities 
in  the  second  set  of  terms,  etc.  An  investigation  of  the  first  order 
singularity  leads  to  its  employment  in  a  superposition  method  to  construct 
various  second  order  singularities. 

Paragraph  5  contains  an  examination  of  the  second  term  of  Eq.  (4.10) 
and  hence,  the  development  of  the  second  order  singularities. 


It  is  shown  in  Paragraph  4  that  the  first  term  of  Eq.  (4.10)  may  be 
thought  of  as  corresponding  to  a  concentrated  load  at  the  origin.  An 
examination  of  the  second  term  in  Eq.  (4.10)  revealed  the  following  facts: 


a)  Each  of  the  diagonal  components  m. .  (no  sum)  produces  a  force 


doublet  along  axis  x. 


li 


b)  Each  of  the  off-diagonal  components  m..(i^jy.lc)  produces  a  force 
couplet  about  the  x^  axis 


c) 


Symmetry  of  oP^ 


(n) 


about  one  or  two  axes  insures  that  certain 


components  of  m„  vanish. 


By  dividing  m. .  into  its  symmetric  and  antisymmetric  parts  m. .  and 


,2S 


ij 


it  is  found  that  the  solution  U7  derived  from  m. .  is  identical  to 

i  x 


m.  . , 

ij 

the  strain  E7^  (no  sum)  corresponding  to  jne  first  order  singularity 

and  represents  a  force  doublet  in  the  x.  direction.  Further,  it  is  noted 
2S  2S  9S  1  1  1  1 

that  +  U,  +  U3  =  +  1^22  +  ^33  anc*  a  center  of  dilatation. 

Next,  it  is  shown  that  the  antisymmetrical  part  of  the  second  order 

singularity  produces  which  is  identical  to  the  rotation  about  the  x. 

1  1  1 

axis  corresponding  to  the  fundamental  singularity  and  represents  a 

force  couplet  about  the  x.  axis.  Finally,  the  Galerkin  vector  producing 

1  3 

a  third  order  singular  solution  is  expressed  as  ..*  j  (r)nu  ^  and  G^(x) 

corresponding  to  a  double  center  of  dilatation  is  presented. 


Next,  an  analogy  is  made  between  the  series  representing  the 
Galerkin  vector  pertaining  to  all  isolated  singularities  and  the  spherical 
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harmonics  series  of  potential  theory.  It  is  found  that  the  G^(x) 

represented  by  '  can  cl°finac*  as  a  spherical  biharmonic. 

Further,  it  is  demonstrated  that,  because  there  are  dependencies  among 

coefficients  of  G^(x)  beginning  with  the  0( •  )  torn,  there  are  an 

infinite  number  of  combinations  which  will  yield  any  given  singular 
solution. 

Finally,  furthering  the  analogy  with  potential  theory,  Betti's 
reciprocal  theorem  is  presented.  Tw»  stress  states  are  chosen,  one 
with  a  zero  body  force  field  and  one  with  an  isolated  singularity  at  the 
origin.  Applying  the  reciprocal  work  theorem,  an  expression  is  presented 
which  enables  one  to  reproduce  any  order  derivatives  of  the  displacement 
field  of  a  regular  stress  state  acting  thr  'ughout  V  if  the  displacements 
and  tractions  are  known  on  V  by  simply  choosing  the  appropriate  isolated 
singular  solution. 

A  physical  interpretation  of  the  most  important  isolated  singu¬ 
larities  (i.c.,  concentrated  force,  force  doublet,  force  couplet, 
center  of  dilatation,  and  double  center  of  dilatation)  has  been  given 
for  isotropy.  The  same  approach  could  undoubtedly  be  used  in  the  case 
of  anisotropy .  As  stated  previously,  proof  that  the  derivatives  of  the 
limit  value  of  the  auxiliary  function  must  be  identical  to  the  limit  of 
the  derivatives  of  the  auxiliary  function  must  he  given. 
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